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Abstract. We study the problem of coexistence in a two-type competition 
model governed by first-passage percolation on TJ^ or on the infinite cluster in 
Bernoulli percolation. Actually, we prove for a large class of ergodic station- 
ary passage times that for distinct points x,y there is a strictly positive 
probability that {z £ U^; d(y, z) < d{x, 2)} and {2 G U^; d{y, z) > d{x, 2)} are 
both infinite sets. We also show that there is a strictly positive probability 
that the graph of time-minimizing path from the origin in first-passage per- 
colation has at least two topological ends. This generalizes results obtained 
by Haggstrom and Pemantle for independent exponential times on the square 
lattice. 



1. Introduction 

The two-type Richardson's model was introduced by Haggstrom and Pemantle 
in [2j as a simple competition model between two infections: on the cubic grid Z'', 
two distinct infections, type 1 and type 2, starting respectively from two distinct 
sources si, S2 € U^, compete to invade the sites of the grid Z'^. Each one progresses 
like a first-passage percolation process on Z'', governed by the same family (t(e))egE<* 
of i.i.d. exponential random variables, indexed by the set E'^ of edges of U^, but 
the two infections interfere in the following way: once a site is infected by the type 
i infection, it remains of type i forever and can not transmit the other infection. 
This leads to two very different possible evolutions of the process: 

• cither one infection surrounds the other one, stops it and then goes on 
infecting the remaining healthy sites as if it was alone, 

• or the two infections grow mutually unboundedly, which is called coexis- 
tence. 

The probability that, given two distinct sources, coexistence occurs is of course 
not full, and the relevant question is to determine whether coexistence occurs with 
positive probability or not. Although this competition problem is interesting in 
its own right, it is also a powerful tool to study the existence of two semi- infinite 
geodesies (or topological ends) of the embedded spanning tree in the related first- 
passage percolation model. Thus, Haggstrom and Pemantle proved that coexistence 
for any two initial sources in the two-type Richardson's model on 1? occurs with 
positive probability, and consequently that in first-passage percolation on 2? with 
i.i.d. exponential passage times on the edges, the probability that there exist at 
least two topological ends in the embedded spanning tree is positive. 
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Their results strongly rely on an interacting particle representation of the prob- 
lem which is typical of the exponential passage times. The aim of this paper is 
to extend these results to more general passage times, where this representation is 
not available anymore or at least much less natural. We consider here stationary 
ergodic first-passage percolation on Z"^', d > 2 (and also on an infinite cluster of 
Bernoulli percolation), and prove that, for any two distinct sources, the probability 
that coexistence occurs is strictly positive. As a consequence, we obtain that in the 
related first-passage percolation on Z"^, the probability that there exist at least two 
topological ends in the embedded spanning tree is positive. 

The structure of the proof is the following. First, the key step is to prove that 
there exist two sources such that coexistence occurs, and this is the aim of Section 3. 
Heuristically, the shape theorem of first-passage percolation, combined with the fact 
that the two infections have the same speed, gives the intuition that the largest the 
distance between the two sources is, the hardest it is for one infection to surround 
the other one. More precisely. Theorem 13 . II savs that if d{x,y) denotes the travel 
time between the sites x and y, then there exists a site x such that the event 

• and the set of sites z such that d(0, z) < d{x, z) is infinite, 

• and the set of sites z such that d(0, z) > d{x, z) is infinite 

has positive probability. The proof of this result relies on the existence of a direc- 
tional asymptotic speed in the related first-passage percolation model. 

The next step is to transfer the coexistence result for these sources to any two 
initial sources; this is done by a modification argument of the configuration around 
the sources using a finite energy property for the passage times. Roughly speak- 
ing, this result expresses the fact that non-coexistence is due to a local advantage 
obtained by one infection at the first moments of the competition. The two topo- 
logical ends result is shown by a similar modification argument. These results are 
proved separately in Section 4 for diffuse passage times and in Section 5 for integer 
passage times. 

The last section is finally devoted to the study of a probabilistic cellular automata 
describing a discrete competition model between two infection types, related to the 
chemical distance in super-critical Bernoulli percolation on W/^ . 

We start now with a reminder of the result of existence of directional asymptotic 
speeds in classical first-passage percolation, and an extension of this result to first- 
passage percolation on an infinite Bernoulli cluster. 

2. Reminder on the directional asymptotic speed results 

In classical first-passage percolation, one has the well-known directional asymp- 
totic speed result: if (t(e))ggEd are i.i.d. non-negative integrable random variables, 
then for every x €1/' , there exists ^{x) > such that a.s : 

t(0,nx) E t{0,nx) , , 
lim ~ lim = [i\x). 

n^QO Ti n—>QO 77, 

This result has been extended in full details in a previous work of the authors |3] 
to first-passage percolation on an infinite Bernoulli cluster. The aim of this section 
is to introduce an adapted framework and to recall, without proofs, the results 
needed in this paper. 
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Grid structure of l/. In the foUowing, d > 2. We denote by the graph whose 
set of vertices is li'^, and where we put a non-oriented edge between each pair {x, y} 
of neighbor points in Z'^, i.e. points whose Euclidean distance is equal to 1. This set 
of edges is denoted by E''. A (simple) path is a sequence 7 = {xi, X2, . ■ ■ ,Xn, Xn+i) 
of distinct points such that Xi and x^+i are neighbors and Ci is the edge between 
Xi and Xi-^i. The number n of edges in 7 is called the length of 7 and is denoted 
by |7l- 

For any set X, and w G Z'', we define the translation operator 9u on by the 
relation 

VW e X Ve e E'' (6l„w)e = OJu.e, 

where u.e denotes the natural action of Z'^ on E'*: if e = {a, 6}, then u.e — {a + 
u,b + u}. 

Assumptions and construction of first-passage percolation. 

Denote by Pc{d) the critical threshold for Bernoulli percolation on the edges E'^ 
of Z'^, and choose p £ {pc, 1]. On ^Ie = {0, 1}'^ , consider the measure P^: 

on Qe = {OAf\ Pp = {pSi + {l-p)Sor^\ 

A point a; in rj^; is a random environment for first-passage percolation. An edge 
e G E'' is said to be open in the environment w if We = 1, and closed otherwise. A 
path is said to be open in the environment lu if all its edges are open in uj. The 
clusters of a environment uj are the connected components of the graph induced 
on Z'* by the open edges in w. As p > Pc{d), there almost surely exists one and 
only one infinite cluster, denoted by Coo- On Qs — {^+)^'' j consider a probability 
measure S^, such that: 

on ils = , §iy is stationary and ergodic 

with respect to the previously introduced family of translations of the grid. We 
suppose moreover that satisfies integrability and dependence conditions: 



(1) m = sup / dS,(77) < 



(2) 



3a > 1, 3A,B >0 such that VA C E"^, [ ^ m > B\A\ ] < 



A 
W 



For instance, if is the product measure , assumption (0) follows from the 
Marcinkiewicz-Zygmund inequality as soon as the passage time of an edge has a 
moment of order strictly greater than 2 - see e.g. Theorem 3.7.8 in [Hj. 

Our probability space will then be = il^; x fig. A point in r2 will be denoted 
{uj,r]), with UJ corresponding to the environment, and r] assigning to each edge a 
non-negative passage time which represents the time needed to cross the edge. The 
final probability is: 

on n = QE x^s, P = Pp ® 
In the context of first-passage percolation, as we are interested in asymptotic 
results concerning travel time from the origin to points that tend to infinity, it is 
natural to condition Pp on the event that is in the infinite cluster: 

Pp(.) =Pp(.lO e Coo) and P = Pp®§^. 
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For B e B{nE), with S C {0 ^ 00} and Pp{B) > 0, we will also define the 
probability measure P_b by 



Examples. The previous assumptions of the generalized first-passage percolation 
model include: 

• The case of classical i.i.d. first-passage percolation: take p = 1, i.e. all the 
edges of Z'' are open, and §^ = , where z/ is a probability measure on 
R+. 

• The case of classical i.i.d. first-passage percolation, but allowing the passage 
times to take the value 00 with positive probability: take pdd) < p < 1, 
a probability measure ly on M^, and set E>i, = v"^^ . This is equivalent 
to consider p ^ 1 and — i''^^ , where i> is a probability measure on 
M-i- U {00} that charges 00 with probability 1 — p. 

• The case of stationary first-passage percolation, as considered by Boivin in 
[T]: take p = 1 and Sj, a stationary probability measure. 

The travel time. The chemical distance D{x,y) between x and y in Z'' only de- 
pends on the Bernoulli percolation structure uj and is defined as follows: D{x, y) [u) — 
inf-y I7I, where the infimum is taken on the set of paths whose extremities are x and 
y and that are open in the environment oj. It is of course only defined when x and 
y are in the same percolation cluster, and represents then the minimal number of 
open edges needed to link x and y in the environment lo. Otherwise, we set by 
convention D{x,y) — +00. 

For {uj,ri) G ^l, and {x,y) e Z"^ x Z'*, we define the travel time from x to y: 

d{x,y){uj,vi) = inf rf(7) = inf Vyye, 
7 7 ^ — ' 

where the infimum is taken on the set of paths whose extremities are x and y 
and that are open in the environment oj. Of course d{x,y) = +00 if and only if 

D{x, y) = +00. 

A path 7 from x to y which realizes the distance d{x, y) is called a finite geodesic. 
An infinite path 7 = (xi)i>o is called a semi-infinite geodesic if (x„, Xn+i, • ■ • , Xp) 
is a finite geodesic for every n < p. 

Directional asymptotic speed results. In classical first-passage percolation, we 
study, for each u G Z'^\{0}, the travel time d{0,nu) as n goes to infinity. Here, 
as all points in Z'^ are not necessarily accessible from 0, we must introduce the 
following definitions: 

Definition 2.1. For each u e Z'^\{0} and B e B{nE), let 

r„^H =inf{n > l;0„„w G B}, 
define the associated random translation operator on fl ~ Qe x 

n-1 

and the composed version T^f ^uj) = ^ T^f ((Of . 

fc=0 
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Note that only depends on the environment w, and not on the passage times 
77, whereas the operator acts on the whole configuration {uj^rf). The next step 
is to study the asymptotic behavior of such quantities: 

Lemma 2.2. 0^ is a -preserving transformation, is ergodicforF and 
E = Tf' = and . 1 ™ - ^ 



Proof. The idea is to prove that classical ergodic theorems can be applied. □ 

We turn now to the study of the quantity analogous to d(0, nu) in the classical 
first-passage percolation: 

Lemma 2.3. Let B e Bins), with B C {0 ^ cx)} andPp{B) > 0. Foru G Z'*\{0}, 
there exists a constant f^ > such that 



n 



The convergence also holds m L1(Pb). Moreover, < E p d{0,T^u) < +00. 

Proof. These results are proved with full details when B = {0 00} in 3 . Since 
the proof is essentially the same, we omit it. □ 

Now, for each u G Z'^\{0}, we define the asymptotic speed in the direction u by 

fi{u) = Pp(0 ^ ^)f^ 
for the choice A = {0 ^ 00}. We also define /i(0) = 0. 

Corollary 2.4. Let B e Bi^E), with B C {0 <^ 00} and Pp{B) > 0. For u e 
Z'*\{0}, we have: 

d{0,{T^Ju;)u)iuj,r,) ^{u) - 

^ z;—, — r Pb a-S. 

n Pp(B) 

d(0,(T„^„Mu)(c^,77) 



IJ,{u) Fb a.s. 



Proof. We use the fact that ( — ^g'"" j , as a subsequence of ( — ^.X'"" 

V / n>0 V "'" / n>0 

admits the same almost sure limit and the lemma □ 

In PI, it has been proved that n enjoys the properties that are usual in classical 
i.i.d. first-passage percolation: /i is a semi-norm. In classical i.i.d. first-passage 
percolation with passage time law it is well-known that /i is a norm as soon as 
i^(0) < Pc(d)- In m, we gave a long discussion about conditions on implying 
that /Lt is a norm. Particularly, if is a product measure v"^"^ , /i is a norm as soon 
as py{Qi) < Pc{d)- 
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3. Coexistence result 

Consider the first-passage percolation model on previously introduced. For 
every pair x and y of distinct points in Z*^, say that the event Coex(a;, y) happens 
if 

{z e d{x, z) < d{y, z)} and {z G Z"^; d{x, z) > d{y, z)} 
are both infinite sets. 

The goal of the paper is to prove that for every pair of distinct points x,y G 
Z"*, P(Coex(x,y)) > 0. Our proofs always require the assumption that ^ is not 
identically null and we guess that this assumption is close to be optimal. Let us 
detail a particular case where /j, = and coexistence never occurs. Suppose that 
d = 2 and E>i, = h'^^ , with pv{0) > Pc{2) = i. In this case, fu, is identically null, as 
previously noted. Consider two distinct points x,y G 1? . Since pi^(O) > Pc(2), there 
almost surely exists an infinite cluster of open edges with passage time zero. It is 
known that in dimension 2, the supercritical infinite cluster almost surely contains 
a circuit that surrounds x and y and disconnects them from infinity - see Harris [H] 
or for instance Grimmett's book Clearly, the points in this circuit are equally 
d-distant from x {resp. y). So, if x reaches the circuit before y, it necessarily also 
reaches every point outside the circuit before y. Similarly, if x and y reach the 
circuit at the same time, all the points outside the circuit will also be reached at 
the same time by x and y. In both cases, coexistence does not occur. 

The next theorem gives conditions that ensures that coexistence possibly occurs 
for some (random) x^y. 

Theorem 3.1. Let d > 2, p > pdd) , E:^ a stationary ergodic probability measure 
on (M+)*' satisfying 0) and 0), and ji be the related semi-norm describing the 
directional asymptotic speeds. 

Let B £ Bins), with S C {0 ^ oo} and fp{B) > 0, and y e Z'K We have: 

z/E d(0, Ti^^y) < then PB(Coex(0, T/'^y)) > 0. 

Moreover, if x £ Z'^ is such that p{x) > 0, then y — rx satisfies to the previous 
condition provided that r is large enough. 

Note that when p = I, which corresponds to classical first-passage percolation, 
we can take B — He, and then T^^y is simply equal to y. 

Before beginning the proof, we want to describe an elementary and clever trick 
used by Pemantle and Haggstrom in 5 that will also be useful here. Consider 
figure O The left picture deals with our problem: if we prove that when M„ goes 

■ • • 

A B O 



Figure 1. The symmetry argument 

to the infinity on the right {resp. on the left), then M„ is infinitely often closer 
{resp. more distant) from B than from A with a probability strictly larger than 
0.5, then coexistence holds with positive probability. 
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Now consider the right picture: for fixed n, {d{0, J„), d{0, /„)) has the same law 
that {d{A, Mn), d{B, Mn j), so the event {d{0, Jn) > d{0, In)} occurs with the same 
probabihty as the event {d{A, M„) > d{B, Mn)}. So, if we show that P(d(0, J„) > 
d(0,/„)) > 1/2 for infinitely many n, the result is proved. 

As Haggstrom and Pemantle said, the idea is that there are sites arbitrarily far 
away from the origin which strongly feel from which source the infection is coming. 
Their modus operandi, in the case of i.i.d. exponentials on Z^, was to control the 
infection rate "from the right to the left" and the infection rate "from the left to the 
right" . The main idea of the proof which follows is that the advantage of the closest 
source can be quantified using the existence of a directional asymptotic speed in 
first-passage percolation. Concretely, we will use the law of d{A, Mn) — d{B, Mn) 
(in fact, the law of d(0, J„) — (i(0,/„)) instead of those of {d{A,Mn) > d{B,Mn)} 
(or {d{0, Jn) > d{0, In)}) to carry the information. 

Proof. {Theorem\EJI Choose y G Z'^\{0} such that 

(3) E<i(0,T=„<|M. 
Let us note 

50 = IW {d{0,z) <d{Ti^^yy,z) < +00} , 

II 2 II 1 — > + oo 

51 = IW {+00 > diO,z) > diTfyy,z)}. 

1 2: II 1 — >+oo 

It is obvious that Coex{0, T^yy) — Sq C] Si. Intuitively, one expects that the 
difference between d(0, z) and d{T^yy, z) will be more important if z G M.y, and 
we will effectively consider such z. For the convenience of the reader, we also note, 
for n e Z+ and x e Z'', fn,x = T^f^x. Define Tq,^ = 0, and for n > 0, 

{Xn = d{0,Tn^y) — d{Ti^y,Tn^y), 
x; = d(fi,„f„,_,)-d(o,f„,_,). 

By triangular inequality, one has |X„| < d{0,Ti^y) and |X,'J < c?(0,ri.j^). 

Note that for oj ^ Sq, X„([x)) < as soon as n is large enough, whereas for 
iu ^ Si, X!^{uj) < for large n. It follows that for uj ^ ^onS"!, X„(cj) -HX4_i(cj) < 
d{0,Ti^y){uj) for large n. Let us define 

Qn = '^i^k + ^fc-i), Zn =^ — and Z = lim 
fe=i " "^+°° 

The previous remark implies easily that 

(4) Vcj^5onS'i Z{uj) <d{0,fi^y){Lu). 

By lemma ITHl c?(0,ri,j^) is integrable under Pb. Since |Z„| < d{0,Ti,y), it follows 
(for instance by Fatou's lemma) that 

E Z = E lim Zn > lim E Zn . 

n — »-|-oo n — >+oo 
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Since d{Ti^y,Tn,y) — d{0,Tn-i.,y) ° 6^ , it follows from the invariance of Pb under 
that 

= ¥, p^d{0,fn^y) -E p^d{0,fn-l^y). 

Then, it follows that E ^^{Xi + X2 + ■ ■ ■ + X„) = E f^d{0, f^^y). Similarly, as 
d{fi^y, fn,^y) = d{0, fn+i.-y) o e^, 

Ep^X = (Ep^d(fi,„f„^_,)-d(0,f„,_,)) 

= E j;^d(0,f'„+i -E j;^d(0,f„-y), 

and E^^{X'q+X[ + -- • + X,;_i) = E j;^d(0, f„,_j,) = E p^d(0, T^^^), using for the 
last equality the fact that Ps is invariant under (O^)" and the fact that a distance 
is symmetric. 

2E p d(0.f„y) __ E j; d{0,fny) 

Then, E p^^n = ^—^^ '—■ Since, via Corollary 12.41 — : — converges 

to ip^^, it follows that 

(5) ErZ>M4- 

Putting together 10), |@J and (0), we see that P_b(S'o n 5i) = - or equivalently 
^b{{So n SiY) = 1 - would yield to a contradiction. This concludes the proof of 
the first assertion. 

The second assertion is a direct consequence of Corollary 12. 41 □ 

One can be a bit perplexed by the fact that the position of the source which may 
coexist with a source at the origin is a random variable. The goal of the next result 
is to come back to deterministic sources. Intuitively, one can guess that the larger 
the distance between the two sources is, the higher the probability of coexistence 
will be. This is the spirit of the next result. 

Theorem 3.2. Under the same assumptions as in Theorem \3.1V suppose moreover 
that /i is not identically null. Then, we have: 

• For X with nix) 7^ 0, there is an infinite set of odd values for n G 
such thatF{Coex{0,nx)) > 0. 

• F{3x,y e 1'^ Coex{x,y)) = 1. 

Let us say a word on the unexpected apparition of odd integers. Of course, the 
result would be the same with the set of integers and generally, this additional 
constraint does not bring much. Nevertheless, one will see later that, in the compe- 
tition context, this additional property sometimes prevents the two infections from 
reaching a point at the very same time; it will also plays a fundamental role in the 
proof of Theorem 15. II 

Proof. Let a; e Z'' be such that nix) > and N e Z_|_. Let A — {0 00} and 
B = An {T^j. is odd}. We have, from the FKG inequalities, 

Fp{B) > Pp{T^^ = 1) = Pp{-x ^ oo)> Pp{-x ^ 00) > 0. 
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By lemma and Theorem 12.41 ^ ^'''^ tends to p^^^, so we can find an odd 

integer r>N with ^ ^'"■'^i'^"") < Mgl. By TheoremO one has P(5'o n^i) > 0. 

By its definition, T^^^^ almost surely takes its values in the set of non-negative 
odd integers. Then, we can write 

k odd 

Then, there exists an odd integer fc e Z+, with Pb(S'o n 5i n {T^,rx = ^1) > 0- So, 
if we note n = kr, we have n> r > N , n is odd and 

P(Coex(0, nx)) > Pp(S)iPi3(5o n 5i n {T]^„ k}) > 0. 

The second point is a consequence of the ergodicity assumption. □ 

4. Mutual unbounded growth and existence of two disjoint 
geodesics for diffuse passage times 

The aim of this section is to prove the possibility of coexistence in general 
two-type first-passage percolation, and to study the existence of two semi-infinite 
geodesies in the corresponding spanning tree. 

We first give our assumptions on the law of the passage times §i/ and define the 
two-type competition model. 

Assumptions. We consider first-passage percolation on Z'', with d>2. The open 
edges are given by a realization of a Bernoulli percolation on the edges E"* of TJ^ 
with parameter p G (pcid), 1]: 

on Qe = {0, lf\ Pp = (pd, + (1 -p)<5o)®^'. 

The passage times of the edges are given by a probability measure Si,: 

on Qs = J is stationary and ergodic. 

Finally, we consider the product measure P = Pp (g) §^ on fi^ x fls- We also need 
two distinct initial sources si,S2 in Z''. 

Lemma 4.1. If x G Z'^ is such that d{si,x) < d{s2,x) then for every y in an 
optimal path realizing d{si,x), we have d{si,y) < d{s2,y). 

Proof. Denote by 7(si,x) an optimal path from si to x, and suppose that there 
exists y G 7(51,0;) such that d{s2,y) < d{si,y). Then, by triangular inequality, 

d{s2,x) < d{s2,y) +d{y,x) < d{si,y) + d{y,x) 

but as ?/ G 7(si, x), d{si, x) = d{si,y) + d{y, x) and then d{s2, x) < d{si,x), which 
is a contradiction. □ 

This allows us to define the following two-type first-passage percolation model. 

Definition 4.2. Under the previous assumptions, we set: 

Ai = {x G Z,'^ , d{si, x) < d{s2,x)}, and A2 — {x G Z,'^ , d{s2, x) < d{si,x)}. 

Ai is the set of sites in that are finally infected by type i infection. The time 
of infection of x G Z'^ is t{x) = ini{d{si,x), 1 < i < 2}. We say that x is finally 
infected if t{x) < 00. 
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Note that the set of finally infected points could he larger than the union of Ai 
and A2 : we can not a priori exclude that a point x could be reached simultaneously 
by the two infections, in which case we call it an infected point without defining an 
infection type. 

We say that the two infections mutually grow unboundedly if the two sets Ai and 
A2 are both infinite. 

Thanks to lemma I^TI the mutual unbounded growth of a two- type first-passage 
percolation starting from si,S2 is equal to the event Coex(si,S2) defined in sec- 
tion 3. Note also that Ai and A2 are connected sets, and if they are infinite, by 
a classical compactness argument, one can find, from each x E Ai, a semi-infinite 
geodesic which is completely in Ai. 

We add in this section assumptions that will ensure the uniqueness of optimal 
paths in first-passage percolation under Pp®§y. If A is a finite subset of E"^, denote 
by ^A<: the (T-algebra generated by {(we, ?7e), e ^ A}. We suppose 

(6) VA finite subset of E'*,Ve G A, Va e R+, §^.(7;^ = a|J^A<=) = 0- 
Lemma 4.3. Under the additional assumption 0), we have: 

i. // 7 and 7' are paths that differ at least from one edge, then 

ii. For every x and y in 'Z'^ , the travel time d{x,y), when it is finite, is a.s. 
realized by a unique path. 

iii. For every a GM., if x,x',y,y' are distinct points in , 
Pid{x,y)-d{x',y')=a)=0. 

iv. // inf d{si,x) is finite, it is realized by a unique source Si. 

Proof. These are classical and not too difficult consequences of assumption . □ 

Now, if t{x) < 00, then x is reached first by a unique infection; the path of 
infection 7(0;) is the unique path from the corresponding source to x that realizes 
t{x). The set of eventually infected points is in this case the union of Ai and A2. 
In other words, we can define uniquely, for each eventually infected point, its type 
of infection and its optimal path. 

The union of (7(a;))j.g2'',t(a;)<oo is then a random forest of two trees T{si) and 
T(s2), respectively rooted at si and S2 and respectively spanning Ai and A2. 

The next theorem ensures the irrelevance of the positions of the two initial sources 
in determining whether mutual unbounded growth occurs with positive probability 
or not. Its proof is based on a modification of the configuration around the sources, 
sufficiently strong to change the initial sources, and sufficiently slight to ensure that 
some geodesies are not modified outside a finite box. 

Lemma 4.4. Consider TH^ , with d > 2 and p G (pc(d),l]. Choose a stationary 
ergodic probability measure S,y on D,s = satisfying to the non-atomic as- 

sumption 1^ and to: 

(7) \/ A finite subset of E'',\fe e A,\/e > 0, S^{rie < e\TA-) > a.s. 
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Ifp — 1, we add the assumption that the support of the passage time is conditionally 
unbounded: 

(8) y A finite subset of E'^,\/e £ A,yM > 0, S^(77e > M|^A<=) > a-^- 
Then if si, S2 and s'l, s'^ are two pairs of distinct points in 11^ , 
P(Coea;(si,S2)) > ^ ^V{Goex{s\, s'^)) > 0. 

Let us comment the two assumptions {TJ and ©. They have the form of finite 
energy properties, which are usual in modification arguments: it enables to force 
the occurrence of a wished event inside a finite box. But they also enable the 
passage time of an edge to take as small - and as large when p = \ - values as 
we like. This is rather a technical assumption that could probably be relaxed. For 
instance, assumptions Q and are satisfied for §^ = v®'^ with V is equivalent 
to Lebesgue's measure on R+. 

The next result says that the mutual unbounded growth in the two-type first- 
passage percolation model and the existence of two distinct semi-infinite geodesies 
in the embedded spanning tree in the corresponding first-passage percolation model 
are equivalent. 

Lemma 4.5. Under the same assumptions as in Theorem \4-.4l 

3si,S2 G Z'' such thatF{Coex{si,S2)) > 
jp / there exist two edge-disjoint semi-infinite geodesies 
\ in the infection tree rooted in 

Combining these results with Theorem 13.21 we obtain: 

Theorem 4.6. Consider , with d > 2 and p £ {pcid)^!]. Choose a stationary 
ergodic probability measure on Vis ^ satisfying to the integrability as- 

sumptions Qj, 0) and such that the related semi-norm ^ describing the directional 
asymptotic speeds is not identically null. 

Suppose moreover that Si, satisfies to the non-atomic assumption 

VA finite subset ofE'^,\/e G A,Va £ M+, §^(r/e = al^A-) = 0, 

and to the finite energy property 

VA finite subset o/E'',Ve G A,Ve > 0, E,„{rie < e\TAc) > a.s. 

Ifp — I, we add the assumption that the support of the passage time is conditionally 
unbounded: 

VA finite subset of¥f\\Je £ A,VAf > 0, S,.(?7e > M\TA_a) > a.s. 

Then, 

1. yxj^ye Z'', P{Coex{x,y)) > 0. 

2 p f there exist two edge-disjoint semi-infinite geodesies 
\ in the infection tree rooted in 

Examples. 1. Consider first-passage percolation on lA, d > 2, with a family 
(i(e))egE<* of i.i.d. non negative random variables with a non-atomic unbounded 
support containing 0, for instance an exponential law as in Richardson's model. 
Then 

• For the two-type competition model, the probability of mutual unbounded 
growth is positive for every pair of distinct sources in Z'*. 
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• For the first-passage percolation model with one source, the probability that 
the embedded spanning tree of has two edge-disjoint infinite branches is 
positive. 

These results were proved by Haggstrom and Pemantle in the pioneer work j2| for 
Richardson's model in dimension 2. Our results positively answer to the questions 
asked by Haggstrom and Pemantle about extensions of their coexistence result to 
higher dimensions and more general distributions for passage times. 

2. Consider first-passage percolation on Z'^, d > 2, with a family (^(e))^^^'' oi 
i.i.d. non negative random variables whose law has no atom excepted in oo {i.e. 
edges can be closed with positive probability) and has in its support, for instance 
t ~ p.W[o,i] + (1 ~ p)-Soo, with Pc{d) < p < 1. Then 

• For the two-type competition model, the probability of mutual unbounded 
growth is positive for every pair of distinct sources in Z''. 

• For the first-passage percolation model with one source, the spanning tree of 
the infinite open cluster has two edge-disjoint infinite branches with positive 
probability. 

Remarks. We evoke here some possible extensions of these results. 

1. In the spirit of Deijfen and Haggstrom's work |2], we could have considered 
competition models with fertile finite sets as initial sources rather than points. As 
the argument is a local modification argument around the sources, our proof can 
be adapted to generalize the irrelevance of the initial sources result: if 81,82 and 
8'i, 82 are two pairs of fertile finite sets in Z'', 

P(Coex(5i, 82)) >0^ P{Coex{8[, 8'2)) > 0. 

2. Let us say a word on multitype first-passage percolation. The definitions 
concerning the two-type first-passage percolation can be generalized in the obvious 
manner to consider a competition model between N infections starting from N 
sources si, S2, sn and trying to invade the sites of . In this context, the event 
Coex(si, S2, Sat) is defined as the event that there finally exist an infinite set 
of infected points of each type of infection. Theorems 14.41 and 14.51 can be proved 
in the same manner for iV-type first-passage percolation. The only difficulty is to 
ensure that the considered initial sources si, S2, sm are susceptible to give rise 
to a coexistence configuration: this means initial sources si, S2, for which it 
is possible to find a family of N infinite paths {Ti)i<i<N such that for every i, Ti 
starts from s,; and such that any two of these paths have no point in common. 

Unfortunately, the coexistence result Theorem 14.61 is not available for N sources, 
as it relies on Theorem 13.11 which is only valid for two sources, and whose proof 
doesn't seem to be easy to adapt to more sources. 

We can now begin the proofs of these results. As the arguments of lemmas [4.41 
and 14.51 are very similar, we give the proof of lemma r4. 41 in full details, and give 
only indications to adapt the proof for the geodesies problem. 

Proof. (lemma \4.4\l - Choose si, S2 and s'^, s'2 two pairs of distinct points in and 
denote by A an hypercubic box in Z'' large enough to contain si, S2 and s'l, S2- We 
also define 9A = {a; e A, 3y ^ A, ||a; - y\\i = 1}. 

By enlarging A if necessary, we can assume that si, S2, s'^, s'2 are at a dis- 
tance at least 3 from 9A. For an edge e £ E'^, we say that e £ A if and only 
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if its two extremities are in A and at least one is not in dA. For a point {uj,r]) in 

n^iiExns = {0,1}^' X (M+f , 

(wa,77a) = {(we,?7e),e e A} and (wa<^ , ?yA<=) = {('^e, f?e), e £ E'^yA}. 

For two point x, y that are in A'^U9A, we define dA" {x, y)(j-^) as the infimum, among 
all the paths 7 from x to y whose edges are not in A, of J2e&j Ve, and, when this 
quantity is finite, 7a <= {x, y) is the only path that realizes this infimum. 
Suppose that P(Coex(si, S2)) > 0. This is equivalent to say that 

T(si) contains a infinite branch Fi starting from si, \ ^ 
T{s2) contains a infinite branch F2 starting from S2- I 

Remember that the box A has been chosen large enough to contain si, S2, so there 
exists on Fi {resp. F2) a last point ri {resp. to be in dA. As dA is finite, there 
must exist two distinct points ri,r2 € dA such that: 



(9) 



^ T{si) contains a infinite branch Fi starting from si \ 
and whose last point in dA is ri, 
T(s2) contains a infinite branch F2 starting from S2 
y and whose last point in dA is 7-2. J 



> 0. 



Now, we must introduce the following events: 
Ci = 



There is a simple path {xj)iyi in A"^ such that \\x\ — ri\\i — 1 
and Vi > 1 , 7a<= (ri , ) = {ri,x\, ...,xl_j^,xl). 

^ _ r There is a simple path {xj)j>i in A'^ such that \\xl — r2||i = 1 
^ ^ 1 and^j >l,-fA.{r2,x]) = {r2,xl,...,x]_^^,x]). 

A] = {dA-ixl^ri) + d{ri,si) < dAc{xl,r2) + d{r2,S2)} , 
A^j = {rfA<=(a;|,''2) + rf(r2,S2) < iiA<=(a;|,ri) + d(ri,si)} . 



Let us prove that © implies: 

(10) P I Ci n C2 n Pi 41 n Pi ^2 1 > 0. 

Indeed, suppose that the event in Q is realized. Then the portion {xj)i>i of Fi 
after ri is a good candidate for Ci. As Fi is a branch of the infection tree T{si), 
one has Vi > 1, 7(^1, xj) = (ri, x\, x}_i, xj), which is stronger than the assertion 
needed for Ci. Next, as ri is in Fi, for every i > 1 we know that ri is a point of 
j{si,xl), and the fact that xj is in T(si) implies that 

d(si,ri) + d{ri,x}) = d{si,x}) < d{s2,x}) < d{s2,r2) + d{r2,x\). 

As ri {resp. r2) is the last point of Fi [resp. F2) to be in A, we have d{ri,x\) = 
dA^iriTx}) {resp. d{r2,x\) = dA<'{r2,xl)), and thus A\ is realized. Doing the 
same for C2 and we see that the event that appears in Q is included in 
Ci n C2 n nj>i A\ n nj>i and thus (O is proved. 
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Figure 2. Modification of the infection trees. On the left, compe- 
tition with two sources S\ (black triangle) and (white triangle): 
the box B is in grey, the squares arc the exiting points of the 
branches from A (black for n and white for r-i) and the circles 
are the visible portions of the infinite branches outside A (black 
for x\, and white for xf, ...,x'^). On the right, the con- 

figuration outside A has not changed, but we forced the branches 
to follow j[ (in black) and 73 (dashed) and thus we changed the 
sources into s[, 



Now, as Ci and C2 are in J^\o , conditioning on J^\c gives: 



Ci n C2 n fl n fl 4 

i>l j>l 



/ dP(wA<=,?7AO-lCi(wA-,?7AO-lc2(wA-,??A<=)-lP i^^^]) £ f] A] H f] A: 

\ i>l .7>1 



where we can also write: 



(11) 



{co,r,)Gf]Alnf]A 



i>\ j>l 

Vi > 1, dA'={xl,ri) -dKa{xl, 7-2) < d{r2,S2) -d{ri,si) 
Vj > 1, rf(r2,S2) -d{ri,si) < rfA<=(a;^,ri) -dAc(a;^,r2) 



J^A" I a.s. 



Define 



mi = mi(wA=,r/A») = sup(dA=(a;- ,ri) - dA»(a;- ,r2)) 

i>l 



m2 = m2{LJAc,'nAc 



inf {dAc {x] , ri ) - dA^ (a;^ , r2 ) ) . 
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Now, we have 

(12) P I Ci n C2 n Pi 41 n Pi n {mi < mz} | > 0. 

Indeed, thanks to (|1HI and to lemma IT^ iii. 

P^^n P A^njmi >m2} Tk- =P(d(r2,S2)-d(ri,si) = mi|J^Ac) = 0a.s. 

and then the probabihties in (|12|l and in (|10(l are equal. It is also easy to see, if 
TOi (uji\<: , 77a<: ) < ^2 (lok" , ?7a^ ) , that we can find ai , 02 , 61 , 62 G such that 

6imi < b2m2, 0,21712 < aimi, bi — ai = 1, and 62 ^ fi2 = 1- 

Define also 

M — M{u!\t:^ri\c) — max {flimi, 627712} 

+ max {dA-'iri, z)} + max {c^ac (r2, 2;)} . 

z£dA z£dA 
dA<:{r\,z)<oo dA<:{r2,z)<co 

Now, we build a set G = G'(aJA<: , '7A<:) of good configurations (waj'/a) inside A, 
depending on the configuration outside A. First, since A has been chosen large 
enough, it is possible to draw with the edges in A, a path 7^ that links s'l to ri 
and a path 72 that links s'2 to r2 such that 7J and 72 have no vertex and no edge 
in common. Denote by |7^| {resp. I72I) the number of edges in j[ {resp. 73). We 
define now G as the set of (cja, '7a) that satisfy the following conditions; 

i. Ve e 7^, = 1 and a2m2/\^'i\ < rje < aimi/\j[\, 

ii. Ve G 72, = 1 and 6iTOi/|7^| < T^e < 62"72/|72L 

iii. • if p < 1, then Ve G A\(7( U 72), lu^ = 0. 
• if p = 1, then Ve G A\(7^ U 73), ?/e > M. 

Under the finite energy assumptions ([TJ and |(SJ, on the event {mi < m2}, we 
have P(G(aJA<:, '7A<:)|-?^A<:) > a.s., so (fl^ imphes 



(13)y lci.lc..l{mi<m.}.F(G|.FAOrfP = P(Ci H C2 H (mi < m2} H G) > 0. 

Let us prove now that on the event Gi n G2 fl {mi < m2} fl G, each of the two 
infection trees T(s'i) and T(s2) contains one infinite branch, or in other words 

Gi n G2 n {mi < m2} n G C Coex(si, S2). 

Suppose then that {uj, ry) G Gi n G2 H {toi < m2} fl G. We have, in the configuration 

• 027112 < d(7i) < flimi thanks to condition i. in the definition of G. 

• 6imi < (^(72) < 627M2 thanks to condition ii. in the definition of G. 

• Thus, by difference and by the choice of ai, 61, a2, 62, we have 
mi = TOi(&i - ai) < ^(72) - d{'-f[) < m2{b2 - 02) = m2- 

Moreover, as soon as a path 7 from s'l to ri differs from 7^ by at least one edge, it 
must use an edge e in A\(7j U72), and this edge is either closed or such that rje > M 
thanks to condition iii. in the definition of G. Thus ^(7) > M > oimi > d{'~f[). 
Then, ^[ is the optimal path from s'l to ri. On the other hand, every path 7 
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from s'2 to ri has to use an edge e in A\(7( U 72), and then by the same argument 
d{'j) > M > aimi > ^(7^), and then d(s2,ri) > d{s'i,ri). Consequently, 

riGT(s'i) and -i[cT{s[) and d{s\,ri) = d{-f[). 

In the same manner, 

r2eT(4) and 7^ C r(4) and ^(4, r2) ^ ^(72). 

Let us prove now that the path (a;J)i>i, given by Ci, is an infinite branch of T(s[). 
Let i > 1, and let 7 be a path from s'l to xj that doesn't exit from A in ri but in 
z ^ r\. If 7 is the minimal path from s'-^ to x\ , we must have 

d{s\,x\) = d(7) = d(si,z) + (iAc(z,a;-) < d{s\,rx) + d[^.(rx,x\) 
and then d{s\^ z) < d{ri, s[) + d\o {ri,xl) — d\<: (z, x]) < d{ri, s'l) + (ri, z). 

The last inequality is just the triangular inequality for dA"- But the path 7 must 
then contain at least one edge e in A\(7^ U72), and so such that rje > M 01 uje — 0, 
and then by definition of M, we must have 

d{s[,z) > M > oiTOi + d\a{ri,z) > d{ri,s[) + dAc(ri, z), 

which contradicts the previous inequality. Thus, the portion of the branch between 
s'l and xj that is in A exits A in ri: it is exactly 7J. It remains to see that the 
part of the branch from s'l to Xi that is not in A is exactly (ri, xi, Xi-i, Xi). But 
this condition is always satisfied in the event Ci. Thus the path {xl)i>i, given by 
Ci, is an infinite branch of T{s[), and in the same manner, the path (a:|)j>i, given 
by C2, is an infinite branch of T{s2). We have thus proved the desired inclusion 
Ci n C2 n {mi < 1122} n G C Coex(si, Sj)- 

Now, H13|l ensures that ¥{Goex{s'i, S2)) > 0, which ends the proof. □ 

Proof. ( Theorem \4-5[ . To prove the direct implication, the proof is the same with 
s'l — s'2 — 0. The only difference is to take the two paths 71,72 rooted both in 0, 
with no other point in common and with no edge in common. 

The converse implication can also be proved by an analogous modification argu- 
ment. □ 



5. Mutual unbounded growth and existence of two distinct 
geodesics for integer passage times 

In the previous section, the law of the passage time of an edge was supposed 
to admit no atom, and thus the minimal paths were unique. However, as seen in 
definition 14.21 it is still possible to define a two- type competition model. Choose 
two distinct sources si, S2 G Z''. We define exactly as previously the sets: 

Ai ^ {x £ Z'^,d{si,x) < d{s2,x)}, and A2 ^ {x e Z'^,d{s2,x) < d{si,x)}, 

and say that coexistence occurs, event denoted by Coex(si,S2), if these two sets 
are infinite. 

As in the preceding section, the geodesies inside Ai and A2 will be useful to 
obtain configurations that allow to prove the irrelevance of the initial sources in 
determining whether coexistence happens with positive probability or not. 
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Assumptions. We consider first-passage percolation on Z"*, with d > 2. The open 
edges are given by a reahzation of a BernoulU percolation on the edges E'* of Z"* 
with parameter p £ {pc{d), 1]: 



The passage times of the edges are given by a probability measure S^, : 
on Qs = (1^+)'^ 7 is stationary and ergodic. 

Finally, we consider the product measure P = Pp ® S,y on Sl^ x fig. We say that 
E>„ satisfies condition {Ha) if 



Note that condition I^Jl was exactly: there exists a > 1 such that {Ha) holds. In 
order to have the coexistence result (Theorem 13. 2|l . we suppose moreover that Si, 
satisfies conditions (QJ and {Ha) for some a > 1. 

Theorem 5.1. Let us suppose that Si, satisfies the previous general assumptions, 
and assume moreover that: 

1. The related semi-norm fi describing the directional asymptotic speeds is not 
identically null. 

2. §y is "discrete" : there exists a subset S o/ Z+ such that §^(5*^ ) = 1. 

3. Su satisfies the following finite energy property: for each finite subset A of E*^ 
and each ca S , we have 



4- If p ~ 1, we add the assumption: S is unbounded. 

5. Some stronger integrability is assumed: suppose that one of the three following 
conditions is fulfilled 

• {Ha) holds for some a > d^ + 2d — 1. 

• p ~ 1 and the passage times of bonds have a moment of order a > d. 

• p — 1, §1, is a product measure and the passage times of bonds have a second 



The last integrability condition is the only one that is specific to the discrete 
case: in the diffuse case of the previous section, we could give to a given edge an 
arbitrary small value thanks to (0 , and there was no need to control the length of an 
optimal path. Here, as passage times are integers, we need a stronger integrability 
assumption that helps to control these paths. The second moment assumption 
is classical in i.i.d. first-passage percolation to ensure the shape theorem -see the 
reference article jZ| lemma 3.5; the {Ha) assumption with a > d is the one used by 
Boivin in ^ for the shape theorem in stationary first-passage percolation. Finally, 
the {Ha) assumption with a > -|- 2d — 1 is the one we use in |2 lemma 3.7 to 
obtain the shape theorem when the edges can be closed. Note that, if is the 
product measure i'®*' , assumption {Ha) follows from the Marcinkiewicz-Zygmund 



on Qe = {0, 1}^ , Pp = {p5i + (1 - p)So) 




S,y{uj\ = caI^aO > S„ a.s. 
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inequality as soon as the passage time of an edge has a moment of order strictly 
greater than 2a. 

In any case, we obtain the following estimate: 

Lemma 5.2. There exists Ki > such that for every a G Z'', we can construct a 
random integer M{a) < +cxd such that 

a oo and y <^ oo and \\y\\i > M{a) d{a,y) < Ki\\y\\i. 

Examples. 

• Take p < 1 and 
As a special case, v 
percolation cluster: 

Corollary 5.3 (Geodesies on a Bernoulli cluster). For eachp > pc, consider 
Bernoulli percolation with parameter p. Then, there almost surely exists a 
point of the infinite cluster from which we can draw two disjoint semi-infinite 
geodesies. 

Proof. The considered event is translation-invariant, so its probability is 
nuU or fuU. By Theorem O with §^ = 6f^\ it can not be null. □ 

• Consider a Poisson point process on MJ^ with an intensity proportional to 
Lebesgue's measure. Let a e Z!J_ and define the passage time r]e by r]e — 
1 + arie, where Ue is the number of obstacles around e, i.e. the number of 
points of the Poisson process which are closer from e than from any other 
edge. 

We can now begin the proof of Theorem 15. II 

Proof. ( Coexistence result). The goal is to prove that for each pair si, S2 of distinct 
sources in Z'', P(Coex(si, S2)) > 0. 

By translation invariance, we can suppose si G Z''\{0} and S2 = 0. Since /x is 
not identically null, we can find x G Z"* such that ||si||i and ||a;||i have the same 
parity and such that ^,{x) ^ 0. Thanks to Theorem 13.21 we can consider an odd 
integer no such that P(Coex(0, noa;)) > 0. Note s'^ = uqx. We are going to prove 
that P(Coex(0,si)) > 0. 

Take Ki > and M{0) and M{s[) as defined in lemma lET^ Since 

lim P({M(0) < n} n {M(s'i) < n} n Coex(0, s[)) = P(Coex(0, s[)) > 0, 

n — >-\-oc 

we can find an integer Ri such that 

P ({A/(0) < Ri} n {M{s[) < Ri} n Coex(0, .s[)) > 0. 

Let A = {a; G Z'^; ||a:||i < R}, for a large integer R whose exact value will be 
fixed later. The idea is then to show that every configuration (w, ry) in the event 
A = {M{0) < Ri} n {M(s'i) < i?i} n Coex(0,si) can be modified inside the baU 
A to get a configuration {lo' ^rj') where Coex(0,si) holds. A classical finite energy 
argument concludes the proof: at first, note that P = Pp ® Si^ also enjoys the finite 
energy property. Now if U is a subset of 17 such that there exists a map f : A ^ B 



— v®^ , where the support oi v \s & finite subset of Z^ . 

— 5i gives the classical chemical distance on a Bernoulli 
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with f{x)A<' = XA<= for each x & A, then P{B) > 0, because 

P(B) = / F{B\J^Ao){x) d¥{x) 
Jn 

> [ F{B\J^Ac){x) dF{x) 

J A 

> I P({/(^)}I^Ac)(x) d^x) 
J A 

> 0. 

Let us explain now the modification inside A. In the following, we will assume 
without loss of generality that the greatest common divisor of the elements of S is 
1. By the lemma of Bezout, we can find a finite family of integers ak and Sk, with 

Sfe G S, such that ^f.akSk = 1- Note 

S+ = {k e S;ak>0}, S- ^ {k G S;ak < 0}, 

bi = smallest odd element of S, 62 = smallest even element of S, 
C = max(Ci, C2), -B = max(6i, 62). 

By convention, if S only contains odd integers, we set 62 = 61 and B = 0. 

The next lemma is a geometrical result, and we omit its proof because it is rather 
tedious and not particularly illuminating: 

Lemma 5.4. Let us consider two fixed points oq, ai G Z'' (not necessarily distinct) 
and two non-negative numbers D and K. Let us note A„ = {x G U^; \\x\\i < n}. 

There exists k = K{ao, ai,D, K) < +00 such that the following holds as soon as 
n > k: 

For each distinct ro,ri G Z'' with ||ro||i — \\ri\\i = n and each integerl which has 
the same parity than \\ao — ai||i and satisfies to \l\ < Kn + D , one can construct 
inside A„ two simple paths 70 from uq to rg and 71 from ai to ri with no common 
point (but maybe ao if oq = ai) and such that: 

\-fo\>Kn + D, \ji\>Kn + D, |7o| - |7i| = /. 

We can now define the radius 

R = max(K(0,Si,B,XiC), Hs'iHi + 2. i?i) 

and define A = Afi. Consider a semi-infinite geodesic starting from {resp. s'l) and 
define by ro {resp. ri) the last point of this semi-infinite geodesic which belongs to 
A. Denote 

L = d(ro,0)-d(n,s[). 

For simplicity, we will suppose, without loss of generality, that L is non-negative. 
Remember that s[ has the same parity than si. Let us define 

J 62 = ^2 and b[ = bi if ||si||i does not have the same parity than L{Ci — C2), 
1^ b[ =b'2 = otherwise. 

and / = L(Ci - C2) + b[ - b'^. 

Note that bi — &2 is odd, unless S only contains odd integers. But in that case, 
Ci — C2 is odd and L has the same parity than ||ro||i + \\ri\\i + \\s[\\i, 
same parity than ||si||i. 
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Thus, I and ||si||i always have the same parity. Note that 

1^1 < \L\\Ci~C2\+B 

< max(d(0, ro),d(si,ri)).max(Ci,C2) + B 

< Ki.m'Ax{\\ro\\i,\\ri\\i).C + B = B + CKiR. 

So, by lemma [5.41 and by the choice we made for R, one can construct inside A 
two simple paths with no common point 70 from to rg, and 71 from si to ri such 
that 

ho\>CKiR + B, I71I >CXii? + S, |7o|-|7i|=^- 

Let us note k — |7o| — {LCi + b[). As proved in the upper bound for |/|, LCi < 
CKiR. Wc thus have 

k > CKiR + B- [LCi + b[) >B-bi>0. 

Obviously, |7o| = LCi + b[ + k and |7i| = LC2 + b'2 + k. Define also the following 
quantity M that will played the role of an "infinite" passage time for open edges: 

M = max aiSi + b[b2 + kbi,L ^ {-ai)si + {b'2 + k)bi 

{ i£S+ leS- 

+ max{d{x,y){OAUJAo,T]), \\x\\i = \\y\\i = R} . 

Note that AI is in T\c, the cr-algebra generated by {(cje, ?7e), e G A"^}. Now define, 
for every {uje^ile) G the configuration (w',?]') € fl: set (1^^,77^) = {uJe,rje) for 
e G E'*\A and define (w^,?/^) inside A as follows: 

i. If p < 1, Ve G A\(7o U 7i), = and 77^ ~ 61, but this value does not play 
a special role; 

if p = 1, Ve G A\(7o U 71), r;^ > M, and uj'^ = 1. 

ii. Ve G 7o U 71, = 1. 

iii. Assign a passage time to edges in 70 as follows (remember that I70I = 
LC\ + b\ + k): first, for each i G S+, give to a^L edges the value t]'^ = Si 
and next complete giving to b'l other edges the value rj'^ = 62 and to k other 
edges the value rj'^ = 61. 

iv. Assign a passage time to edges in 71 as follows (remember that |7i| = 
LC2 + 62 + first, for each i £ give to —atL edges the value r]'^ = 
and next complete giving to the remaining b'2 + k edges the value rj'^ — bi. 

Now we immediately obtain: 

v'e = L ajSj + b'ib2 + kbi, 

ee7o ies+ 

ee7i ieS- 

Y.'^'e-Y.'^'e = i(^a,s,) +^i&2-?>2fei =i = d(0,ro)-d(s'i,ri). 

ee7o eG7i ieS 

For (WjTy) G A C Coex(0,S2), there exist two infinite geodesies starting from 
and s'-y. Let us denote by Fq {resp. Fi) the part beginning at tq [resp. ri) in the 
geodesic starting from {resp. s'l) in the configuration [uj^if). 
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We are going to prove that 70 UFq {resp. 71 UFi) is an infinite geodesic starting 
from (resp. si) in the configuration (to' , 77'). Let x be a point of Tq. Let us prove 
that an optimal path from to a; in the configuration (w', 77') is included in 70 UFq. 

Let 7 be an optimal path from to x in the configuration [uj'^rj'), and denote 
by z the point from which the path 7 exits from A. We have 

But since ?/(e) = i 0,^5^+6']^ &2+^&i and c?(ro, 77') < (i(ro, 2:)(0Aa;Ac, ry), 

eG7o iSS+ 

it follows that (i(0, z) < A/. By definition of Af , it ensures that 7 does not use any 
bond in A, except those used in 70 U 71, and particularly, it implies that z = rp, 
and thus an optimal path from to x is included in 70 U Fq. 

Similarly, let 7 be an optimal path from si to x in the configuration [uj' , rj'), and 
denote by z the point from which the path 7 exits from A. We have 

eS7 eG7i 

But since X;ee7i '7'(e) = L X^^eS- (^"^)*»+(^2+^)^i and (i(ri , z)(cj', 77') < (i(ri, z)(Oalja 
it follows that d(si, z) < M. By definition of M, it ensures that 7 do not use any 
bond in A, except those used in 70 U 71, and particularly, it implies that z = ri, 
and thus an optimal path from si to x uses 71 to exit from A. 
Let us now prove that if a; € Fq, d(0, x){uj' , rj') < c?(si, x){uj' , rj'). 

d{si,x){u}' ,ri') = d(si,ri)(u;',7/) + d(ri,a;)(0AWA=,?7')i 

= d{si,ri){uj' ,1]') + d(ri,x)(OAWA<=,?y), 

d(0,x)(w',r;') = d(0,ro)(c^',r;') + d(ro,a:)(OAC^A=,77'), 

= d(0,ro)(w',?7') + d(ro,a;)(OAWA<=,?y)- 

Consequently, 

(d(si,x)-d(0,x))(c^',7/) 

= ((i(si,ri) - (i(0,ro))(a;',77') + (i(ri, a;)(OAWA<^, ??) - d{ra,x){Qj^u}K^ ,vi), 
= {^{s^.ri) - d{0,ro)){uj,r]) + d(ri , x) (Oawa<= , r?) - d(ro, a;)(OAt^A<=, 
> d{s[,x){u;,r]) -diO,x){uj,r]) >0, 

because a; G Fq, which is a part of the infinite geodesic issued from in the configu- 
ration (w, rj). Thus 7o UFq is an infinite geodesic issued from in the configuration 

In the same manner, working symmetrically with Fi, we prove that 71 UFi is an 
infinite geodesic issued from si in the configuration (a;', 77'), and finally {uj',ri') G 
Coex(0,si). □ 

Proof. (Geodesies result). The goal here is to prove that 

p / there exists two distinct semi-infinite geodesies 
I starting from for the random distance d 

The proof is exactly the same as the previous one. The only difference is to use the 
single source rather than two distinct sources 0,s[. The geometrical structure 
of the modification is once again given by lemma 15.41 and the adjustment of the 
values is made as before. □ 
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As seen previously, a trouble with integer passage times is that some points can 
be reached at the very same moment by the two distinct infections. This case can 
be ruled out under some extra assumptions, and this is the goal of the next result. 
But first, for two distinct sources x,y Q Z'', we say that the event Sep-Coex(a;, y) 
happens if 

{z G Z'', d{x, z) < d{y, z)} is infinite 
and {z G Z'', d{x, z) > d{y, z)} is infinite 
and Vz G Z'^, d{x, z) ^ d{y, z) unless d{x, z) = d{y, z) — +00. 

We have the following result: 

Lemma 5.5. Denote by O the set of non-negative odd integers, and as previously, 
let d > 2, p > Pc{d), §1/ a stationary ergodic probability measure on satisfying 
0J) and 0). Then, for a; G Z'' with \\x\\i odd, 

¥{Coex{0, x)\Sep-Coex{0, x)) = 0. 

Proof. By the assumption we made on /i, the length of a path from x to y has the 
same parity than ||a; — y\\i. So, the identity d{x,z) = d{y,z) can only happens if 
\\x — 2/11 1 is even. □ 

Now, for a given point x with ||a;||i odd, the fact that P(Scp-Coex(0, a;)) > can 
be obtained as a consequence of Tlieorem l3.2l or Theorem 15. II Note also that when 
the assumptions of lemma are fulfilled, is always a norm: since the passage 
time of a bond is an odd integer, it is at least equal to 1. Then, it is easy to see 
that for each x G Z'', we have [i(x) > \\x\\i. 

6. An EXAMPLE OF A DISCRETE TIME COMPETING PROCESS 

The last section is devoted to the study of a natural example of a non-trivial 
dynamical system which can be studied with the help of Theorem l5.1l and lemma l^31 

Consider two species, say blue and yellow, which attempt to conquer the space 
Z"*. At each instant, each fertile cell tries to contaminate each of its non-occupied 
neighbors. It succeeds with probability p. In case of success, the non-occupied cell 
takes the color of the infector. If a yellow cell and blue cell simultaneously succeed 
in contaminating a given cell, this one takes the green color. If a green cell and 
another cell simultaneously succeed in contaminating a given cell, this one takes 
the green color. At the next step, the individuals that have just been generated 
are fertile, but the previous generation is no more fertile. We make the following 
assumptions: 

• the success of each attempt of contamination at a given time does not 
depend on the past, 

• the successes of simultaneous attempts to contamination are independent. 
The first assumption allow a modelization by an homogeneous Markov chain. 

Markov chains satisfying to the second condition are sometimes called Probabilistic 
Cellular Automata (PC A). 
Let us define 

S = {0, blue, yellow, green, blue* , yellow* , green* } , 

where is the state of an empty cell, blue, yellow, green the states of fertile cells, 
and blue* , yellow* , green* the states of unfertile cells. 
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Since we will study the evolution of a system which starts with only two cells, we 
will only deal with configurations in which a finite numbers of cells are nonempty. 
So, we will deal with a classical Markov chain on the denumerable set 

C = {(, & S'^" ; 3A finite, (,k = Q iov k eZ^K}. 

We now define for color G Act ~ {blue, yellow, green} : 

n(color,a:)(0 = \{y € Z'*; ||a; - y\\\ = 1 and = color}], 

and 

s(color,a;) = 1 - (1 

which represents the probability that at least one neighbor of x succeeds in infecting 
X with the given color. 

The considered dynamics form an homogeneous PCA with space state S and 
whose local evolution rules are given by 

^2 _ pyi(ycllow,x)+n(blue,a;)+n(green,a;) J£ ^ _ ^ _ Q 

s(blue, x){l - p)»(yeiiow,x)+n(green,x) if s = and t = blue 

s(yellow, x){l - p)"(biuc,:^)+n(grocn,i.) jf ^ ^ Q ^nd i = yellow 

Px{s, t) = \ s(green, a;) + (1 — s(green, a;))s(blue, a;)s(yellow, x) if s = and t = green 
1 if s S {blue, yellow, green} and t = s* 

1 if s e {blue*, yellow*, green*} and t = s 

otherwise. 

In term of Markov chains, it means that the transition matrix is defined by 

The product is convergent because only a finite numbers of terms difii'ers from 1. 
With the help of the tools that we have developed above, we will prove the 

following theorem: 

Theorem 6.1. Letp > pc- For SypUow, Sbiue G '^'^ with s^iue Syeiiow, let us denote 
by ]Pj),sjeiiou„s(,i„r the law of a PCA {Xn)n>o following the dynamics described above, 
and starting a confi,guration with exactly two non-empty cells: a blue cell at site 
Sblue, 0, yellow cell at site Syeiiow, the others cells being empty. Then, 

Pp,s,e»„»,s«.e(Vn e Z+ 3{x,y) e Z'^ x Z^, X„(a;) = blue and X„(y) = yellow) > 0. 

If moreover \\syeUow — SbiueWi is odd, green cells never appear. 

The following lemma gives the link between this PCA and our competing model. 

Lemma 6.2. Consider a probability space where lives a family (a;e)eeE'' of inde- 
pendent Bernoulli variables with parameter p, which defines a random chemical 
distance D. 
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Figure 3. Two samples of simulation of the competing process 
when p = 0.6. The process is stopped when the border of a given 
box is attained by one of the two species. The color in the picture 
is determined by the time of coloring and the type of the cell. 



Let Syeiiow, Sbiue £ With s blue 7^ SyMow Define 



( blue if D{sbiue, x) < D{syeiiow, x) 

yellow if D{syeiiow, x) < D{sbiue, x) 

green if n ^ D{syeiiow, x) = D{sbiue, x) 

Xn{x) = {blue* if D{sbiue,x) < min{D{syeaow,x),n) 

yellow* if D{syeUow, x) < miii{D{sbiue, x), n) 

green* if D{syeUow, x) = D{sbiue, x) < n 

otherwise. 

Then, {Xn)n>o is an homogeneous PCA with space state 

S — {0, Wite, yellow, green, blue* , yellow* , green*} 

associated to the probabilities Px{s,t) defined above. 

Proof. Let us consider the map 

^ {fx{(.x,(^))xez'i 

where fx ■ S x He ^ S is defined by 

fx{s,uj) — s* for each s e {blue, yellow, green} 



fx{s,u)) 



s for each s e {blue* , yellow* , green* } 



fxiO,t) = { 



blue if Act n {C 

yellow if Act n {^y; ||a; 

green if Act n {^j^; II a; 

green if Act n {^j^; ||x 

otherwise. 



y\ 


1 


= 1 and 


^{x,y} 


y\ 


1 


= 1 and 


^{x,y} 


y\ 


1 


= 1 and 


^{x,y} 


y\ 


1 


= 1 and 


'^{x,y} 



1} = {blue} 
1} = {yellow} 
1} = {blue, yellow} 
1} D {green} 
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By considering Dijkstra's algorithm in the particular case where the travel times 
are constant, it is not difficult to see that {Xn)n>o satisfy to the recurrence formula 
Xn+i = f{Xn,uj). To recognize (^Ti)n>o as a convenient PCA, we will build a 
coupling of Lo with a i.i.d. sequence (w")n>i to obtain the canonical Markov Chain 
representation Xn+i = f{Xn,uj^^)- 

Let (r2, !F, P) be a probability space with a;°, lo^jUj"^, . . . independent {0, 1}^ 
valued variables with Ber(p)®'^ as common law. 

We define Aq — {sbiuo, Sydiow} and recursively 

{Bn+i= {yeZVn 3x^dAr,:\\x-y\\i = lanALo'l^y^ = l} 

[ Ai+l = An U Bn+l- 

Note that the random set Bn^i is measurable with respect to the cr-algebra gener- 
ated by w^, . . . , w"). We define C" recursively by 



I w^+i if e = {x, y} with (x, y) G dA^ x U^XA^ 
I C" otherwise. 



By natural induction, we prove that the law of under P is Ber(p)®'^ . By 
construction, each bond e writes e = {x, y} with (a;, y) E dAn x Z'*\A„ for at most 
one value of n. It follows that the sequence converges in the product topology. 
Let us denote by uj°° its limit. Since the law of C" under P is Ber(p)'*'^ , it follows 
that the law of ijJ°° under P is also Ber(p)'*'^''. 

Now, it is not difficult to see that sequence (X„)„>o defined from uj"^ as previ- 
ously, satisfies to the recurrence formula Xn+i = f{Xn,oj°°), but also to Xn+i = 

It is now proved that {Xn)n>a is an homogeneous Markov chain. The recognition 
of the transition matrix follows from an elementary calculus. □ 

We can now prove the theorem announced above. 

Proof. Clearly, lemma \^7B connects the considered PCA with the random distance 
studied in Theorem 15.11 Here, the passage times of open bonds are identically 
equal to 1, which is obviously an odd number. By lemma this prevents from 
the appearance of green cells when ||syoUow — Sbiuc||i is odd. □ 

Remarks. If ||syciiow — Sbiuc||i ^ is even and if the two species infinitely grow, 
there are necessarily green cells at the boundary between blue cells and yellow cells. 

A natural question is the following: is it possible to have an infinite set of green 
cells surrounding the blue cells and the yellow cells ? The answer is yes, as soon as 

ll^yellow -^bluc 

111 7^ is even: consider figure El 

The picture describes a particular case when = 2, but the reasoning can 
obviously be generalized. 

In this case, the yellow fiow and the blue flow immediately converge to engender 
a green fiow. They also do not develop themselves elsewhere. If the point labelled 
belongs to the infinite cluster, then the result is proved. 

It is now to see that, conditionally to the states of the bonds imposed by this 
picture, the probability that belongs to the infinite cluster is strictly positive, 
which follows from a classical modification argument. 
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Figure 4. Green surrounding the two sources 
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